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Abstract 

We study the degrees of freedom (DOF) of a multiple-input multiple-output (MIMO) two-way X 
relay channel, where there are two groups of source nodes and one relay node, each equipped with 
multiple antennas, and each of the two source nodes in one group exchanges independent messages 
with the two source nodes in the other group via the relay node. It is assumed that every source node is 
equipped with M antennas while the relay is equipped with N antennas. We first show that the upper 
bound on the total DOF for this network is 2 min {2M, N} and then focus on the case of N < 2M so 
that the DOF is upper bounded by the number of antennas at the relay. By applying signal alignment 
for network coding and joint transceiver design for interference cancellation, we show that this upper 
bound can be achieved when N < L^f^J- We also show that with signal alignment only but no joint 
transceiver design, the upper bound is achievable when N < L^J- Simulation results are provided to 
corroborate the theoretical results and to demonstrate the performance of the proposed scheme in the 
finite signal-to-noise ratio regime. 

Index Terms 

MIMO X channel, relay, two-way communication, signal alignment, joint transceiver design. 

I. Introduction 

Wireless communication has been advancing at an exponential rate, propelled by the ever- 
increasing demands for wireless multimedia services. This, in turn, necessitates the development 

Z. Xiang, M. Tao and J. Mo are with the Department of Electronic Engineering at Shanghai Jiao Tong University, Shanghai, 
200240, P. R. China. Email: {7222838, mjh, mxtao}@ sjtu.edu.cn. 

X. Wang is with the Department of Electrical Engineering at Columbia University, New York, USA. Email: 
wangx@ee.columbia.edu. 



2 

of novel signaling techniques with high spectrum efficiency and capacity. Among those factors 
limiting the capacity of wireless networks, interference has been considered as a key bottleneck. 
Recently, two advanced signaling schemes have been proposed to cope with interference and to 
enhance spectrum efficiency: network coding and interference alignment. 

Network coding was originally proposed in 0]| to achieve the max-flow bound for the wireline 
network. The key idea of network coding is to let an intermediate node combine the messages 
it receives and forward the mixture to several destinations simultaneously. Compared with the 
conventional time-sharing based schemes where different destinations are served at different 
time slots, network coding can increase the overall throughput significantly. The first wireless 
application of network coding was the two-way relay channel, where two source nodes exchange 
information with the help of a relay (sometimes referred to as physical layer network coding) 
[El 0. By applying physical layer network coding at the relay, the spectrum efficiency of the 
two-way relay channel can be doubled compared with the conventional schemes. Physical layer 
network coding has also been applied to several other relay-aided wireless networks such as 
multiuser two-way relay networks [|4]|- 10, multipair two-way relay channels 0- flU and multi- 
way relay networks [ITOll - 021. 

Interference alignment was first proposed in lfT3l lfl4l to achieve the maximum degrees of 
freedom (DOF) for the multiple-input multiple-output (MIMO) X channel. It has been shown 
that for the MIMO X channel with every node equipped with M antennas, its total DOF is ^p. 
The key idea is to align the interference signals so that they occupy the smallest signal space, 
leaving more free space for the useful signals. It was shown in 0"5l that the capacity of a iT-user 
time-varying interference channel is characterized by 

C(SNR) = — log(SNR) + o(log(SNR)). (1) 

Thus, independent of the network size, it is theoretically possible that each user achieves half 
the DOF of an interference-free system. Hence interference is not a fundamental limitation 
for such networks. A number of interference alignment schemes have been proposed, such as 
distributed interference alignment, ergodic alignment and blind interference alignment ifToTl - [fT8l . 
An overview on various interference alignment techniques is given in |fT9l . 

Based on the concept of interference alignment, signal alignment was proposed in ll20l to 
solve the network information flow problem for the MIMO Y channel, where there are three 
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users and a single relay, and each user sends information to the other two users via the relay. 
Unlike interference alignment, the goal of signal alignment is to align the signal streams for 
different user pairs at the relay. Combined with network coding, it can significantly increase the 
network's throughput. In 11211 11221 . signal alignment was applied to the generalized K-user Y 
channel. 

In this paper, we consider the network information flow problem for the MIMO two-way X 
relay channel and analyze its total DOR In this network, there are two groups of source nodes 
with each group consisting of two nodes, and a relay node. Each source node in one group 
exchanges independent messages with the two source nodes in the other group with the help of 
the common relay. It is assumed that every source node is equipped with M antennas and the relay 
node is equipped with N antennas. We first show that the DOF of this network is upper bounded 
by 2min{2M, N}. By combining the techniques of signal alignment for network coding and 
joint transceiver design for interference cancellation, we then propose an efficient transmission 
scheme and show that this scheme achieves the upper bound when N < |_^irj • We also show 
that with signal alignment only but no joint transceiver design, the upper bound is achievable 
when N < [_^rJ • Note that the MIMO two-way X relay channel has been considered in ||23l for 
a special case of M = 3, N = 4. In this paper, we consider the general case with arbitrary M 
and N. Moreover, when N > our proposed scheme outperforms the generalized version 

of the scheme in 11231 . 

The remainder of this paper is organized as follows. In Section II, the system model of 
the MIMO two-way X relay channel is described. In Section III, we derive an upper bound 
on the DOF of this channel. In Section IV, we present an efficient transmission scheme and 
give a necessary condition, i.e., iV < [^J , for this scheme to achieve the upper bound. In 
Section V, we show that the necessary condition is also sufficient. In Section VI, we consider a 
special variate of our proposed transmission scheme which reduces to the method in ||23l when 
M = 3, N = 4. Simulation results are provided in Section VII. Finally, Section VIII concludes 
the paper. 

Notations: Boldface uppercase letters denote matrices and boldface lowercase letters denote 
vectors. R, C and Z + denote the sets of real numbers, complex numbers, and positive integers, 
respectively. = max{n e Z + |n < x}. (-) T , (-) H , and Tr{-} are the transpose, Hermitian 
transpose, Moore-Penrose pseudoinverse and trace operators, respectively. E(-) is the expectation 
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operator. Span(H) and Null(H) stand for the column space and the null space of the matrix H, 
respectively. dim(H) denotes the dimension of the column space of H. I N denotes the N x N 
identity matrix and © is the exclusive-OR operator. 

II. MIMO Two-Way X Relay Channel 

Consider a MIMO two-way X relay channel shown in Fig. Q] The channel consists of four 
source nodes with M antennas each and a relay with N antennas. Each source node i, for 
i = 1,2 on the left-hand side (LHS) needs to send an independent message, denoted as 
to each source node i', for i! = 3,4 on the right-hand side (RHS) via the relay. So does each 
source node on the RHS. 

The transmission is implemented in two phases. In the multiple-access (MAC) phase, all four 
source nodes transmit their signals to the relay. The received signal at the relay is given by 

4 

y r = ^2 Hi, r Xj + n r (2) 

i=l 

where y r and n r denote the N x 1 received signal vector and the additive white Gaussian noise 
(AWGN) vector at the relay, respectively; Xj is the Mxl transmitted signal vector by source 
node i with the power constraint E (Tr{xjX^}) < P; H i;r is the N x M channel matrix from 
source node i to the relay. The entries of the channel matrices H i r for, % = 1, 2, 3, 4, and those 
of the noise vector n r , are independent and identically distributed (i.i.d.) zero-mean complex 
Gaussian random variables with unit variance, i.e., CJ\f(Q, 1). Hence, all channel matrices are of 
full rank with probability 1. 

After receiving the signals from the source nodes, the relay forms a new signal x r . and 
broadcasts it to all source nodes, which is known as the broadcast (BC) phase. The received 
signal at the ith source node is given by 

Yi = H n iX r + n i; i = 1, 2, 3, 4 (3) 

where and rij denote the Mxl received signal vector and the AWGN vector at the ith 
source node, respectively; x r is the iV x 1 transmitted signal vector by the relay with the power 
constraint E (Tr{x r x^}) < P; H r i is the M x N channel matrix from the relay to source 
node i. Similar to the MAC phase, we assume that H r j and rij contain i.i.d. CAf(0, 1) random 
variables. 
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Throughout this paper, it is assumed that perfect channel state information (CSI) is available 
at all source nodes and the relay. Additionally, we assume that the source nodes and the relay 
operate in full-duplex mode. 

We define the total DOF of the above MIMO two-way relay X channel as 

d = dx j3 + di A + a? 2 ,3 + d 2) 4 + d 3jl + d 3)2 + d*,i + (4,2 
#(SNR) 

= ] im — ^ — (4) 

SNR^oolog(SNR) 

where dij is the DOF from source node i to source node j, and R (SNR) is the sum rate as a 
function of SNR, where SNR is defined as SNR = P since the noise samples are assumed to 
have unit variance. 



III. An Upper Bound on DOF 

In this section, we derive an upper bound on the DOF of the MIMO two-way X relay channel. 

Theorem 1: Consider a MIMO two-way X relay channel with M antennas at every source 
node and N antennas at the relay. The total number of DOF is upper bounded by 2 min {2M, N}, 
i.e., 

d < 2 min {2M, N} . (5) 

Proof: We first consider the network information flow of one direction, i.e., from source 
nodes 1,2 to source nodes 3,4 via the relay. Then the signal model in (O for the MAC phase 
can be written as 

y' r = Hi jr xi + H 2 ,r.x 2 + n' r 

+n' r . (6) 

Xl2 

Similarly, the BC phase in © can be written as 

y ; = H r X + < * = 3 ,4 (7) 

Applying the cut-set theorem [|24l to each phase, as shown in Fig. [2l the cut-set bound on the 
rate is then given by 

Ri,s + Ria + #2,3 + #2,4 < min {/ (y; ; x l9 x 2 ) , / (y' 3 , y' 4 ; x^,)} (8) 




Xl 

x 2 
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where R it j denotes the information rate from source node i to source nodes j; and I(y; x) 
denotes the mutual information between x and y. For the first term on the RHS of ([8]), we can 
upper bound the mutual information by considering source nodes 1 and 2 as fully cooperating, 
i.e., 

SNR 



/(y^XijXa) < log 



det Ijv + 



min {2M, N} 



Hl2,rH 12 , 



min{2A/,7V} 
i=l 



SNR - A 



12,r 



(9) 



min {2M, N} 

where \] 2 ' r denotes the ith eigenvalue of the matrix H 12)r Hf| r . For the second term on the RHS 
of ©, we first rewrite © as 

y 3 4 = H ri3 4x' r + n' 34 (10) 



where y 34 
have 



/ A 





T 


/ T ,T~ 


/ T ,T 

y 3 >y 4 


1 A 
» n 34 = 


n 3 > n 4 



and H r)34 = [H^ 3 ,H^ 4 ] T . Then similarly we 



min{2M,iV} 



/(y 3) y 4 ;x;)< £ log i + - 



i=l 



SNR ■ AI 



r,34 



min {2M, N} 



(ID 



where A[' 34 denotes the zth eigenvalue of the matrix H r34 H^3 4 . Based on the above results, we 
have 



^1,3 + "1,4 + ^2,3 + ^2,4 

#1,3 + #1,4 + #2,3 + #2,4 



lim 

SNR^oo 



log (SNR) 



- • } v J (yr; x i 5 x 2) ,. ^(y3>y 4 ; x ^ 

< mm < hm — — — — — — , hm , :£ 7 — ^— 

SNR^oo log(SNR) 'SNR^oo log(SNR) 



For the first limit in (fl~2l) . we have 

J(y£;xi,x 2 ) 



lim 

SNR^oo log(SNR) 



< lim 

SNR->oo 



lim 

SNR-s>oo 



min{2M,A/"} 

E log ( i + 

i=i 



SNR-A.^' r 
min{2M,Af} 



log (SNR) 

min{2M,7V} f 

E log (SNR) + log 

1=1 L 



min{2Af,AT} 



log (SNR) 



min {2M, . 



(12) 



(13) 
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For the second limit in (fT2l) . similarly we have 




< min {2M, N} . 



(14) 



Hence we obtain 



di, 3 + di,4 + d 2 ,s + d 2A < min {2M, . 



(15) 



Considering the other direction of the network information flow, we can similarly obtain 



From the above result, we can see that when iV < 2M, the total DOF for the MIMO two-way 
X relay channel is upper bounded by the number of antennas at the relay, which is therefore the 
bottleneck for the spectrum efficiency of the network. In the remainder of the paper, we assume 
that iV < 2M so that the upper bound on the DOF is 2N. Since the transmission scheme for 
the case of iV > 2M will be completely different from that for the case of N < 2M, we will 
leave the case of iV > 2M to future work. 



For relay-aided bidirectional channels, applying physical layer network coding at the relay can 
significantly improve the system's spectrum efficiency; and for multiuser channels, beamforming 
is typically employed to nulled out the multiuser interference. In this section, by applying 
signal alignment for network coding and joint transceiver design for interference cancellation, 
we propose a novel transmission scheme, named as "Signal Alignment with Joint Interference 
Cancellation (SAJIC)" for the MIMO two-way relay X channel to maximize its total DOF. 

A. A Motivating Example for M = 5, N = 8 

As an example, we consider a system where each source node has M = 5 antennas and 
the relay has N = 8 antennas. For this system, the proposed transmission scheme achieves 
di,3 = g?i,4 = ^2,3 = ^2,4 = o?3,i = d 3 2 = c?4,i = ^4,2 = 2. In particular, source node 1 transmits 
codewords s\ 3 , 3 (s\ 4 , sf 4 ) for message Wi i3 (W 1:4 ) by using beamforming vectors v\ 3 , 3 
(vj 4 , Vj 4 ), respectively to source node 3 (source node 4) via the relay. Similarly for the other 
three source nodes. 



d 3 ,i + 4,2 + d 4>1 + d 4 , 2 < min {2M, N} . 



(16) 



Combining ©, (fl~5l) and (fl~6l) . we conclude © which completes the proof. 



IV. Efficient Transmission Scheme 



s 



Step 1: Signal alignment during the MAC phase 

During the MAC phase, there are totally 16 data streams arriving at the relay. Since the relay 
has only 8 antennas, it is impossible for it to decode all the 16 data streams. However, based on 
the idea of physical layer network coding, the relay node only needs to decode some mixtures of 
the symbols. Specifically the key point of the proposed scheme is to obtain the network coded 
messages W 1)3 © W 3>1 , W 1A © W 4 ,i, W 2 , 3 © W 3 , 2 and W 2 ,a © W 4 , 2 at the relay (Note that each 
message consists of two streams.). Inspired by the signal alignment for network coding 11201 . we 
design the beamformers so that the two desired signals for network coding are aligned within 
the same spatial dimension. Taking source node 1 as an example, we align its transmitted data 
streams with the streams from source node 3,4 as follows 

span (Hi ir .v{ 3 ) = span (H^v^) = g* 
span (H 1)r v? 3 ) = span (H^v^) = 
span (Hi^vJ 4 ) = span (H 4 ,rV 4jl ) = g! 

span (Hi, r v? j4 ) = span (l^rvf^) = gj (17) 

where g*, g^, gf, g^ are the signal vectors seen by the relay. Fig. [3] illustrates the notion of the 
signal alignment in the MAC phase where it is seen that there are 8 network coded symbols 
aligned along 8 signal vectors, respectively. With AT = 8 antennas, the relay can then obtain the 
above 8 network coded symbols. 
Step 2: Joint transceiver design for interference cancellation during the BC phase 
During the BC phase, the relay broadcasts these four network coded messages using beam- 
formers u*, u®. More specifically, [uj, u^], [u£, u£], [u^, u£], [uj, u®] are for messages Wi >3 © 
^3,1, Wi,4 © W4,i, W2,3 © ^3,2, 1^2,4 © ^4,2 , respectively. Note that at the receiver side each 
source node suffers from two sources of interference and each component of the transmitted 
signal causes interference to two source nodes. For instance, source node 1 suffers from the 
interference caused by [uj?,u^], [u£, u®]; and [uj, uj!] causes interference to source nodes 2,4. 
Each source node employs a linear receiver, with the 4x5 receiving filter matrix denoted as 
Dj, for i = 1,2,3,4. Denote the effective channel matrix from the relay to source node i as 
H r j = DjH r j. The goal of transmit beamformer design at the relay is to make each component 
of the transmitted signal to lie in the null space of the effective channel matrices of those 
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unintended source nodes. For example, [uj, u 2 ] should satisfy the following condition 



/ 



[ul,u 2 r ] C Null 



H rj2 

H r 4 



Null 



d2 )3 H ri 2 
d2 )3 H ri 2 

^2,4^,2 
d 4,l H r,4 
d 4 2 Hr,4 



\ L ^4,2^,4 



(18) 



/ 



where d^ ■ is the lxM receiving filter vector for source node % to extract the A;th data stream 
from source node j. 

However, the dimension of the matrix on the RHS of (fT8~l) is 8 x 8 and in general it is full- 
rank if each source node chooses its receiving vectors {d^ •} independently. Therefore it is not 
possible to find [u} r , u 2 ] that satisfies (fT8~l) . In order to circumvent this problem, we consider joint 
transceiver design for the source nodes and the relay. More specifically, for source nodes 2,4, 
we constrain their receiving vectors to satisfy 



span ((d^ >4 H ri2 ) T ) = span fjd^H^)' 
span ((d2 i4 H ri2 ) T ) = span ((d^H^)' 



(19) 



That is, the effective channel matrices for the network coded message Wz^ © from the 
relay to source nodes 2 and 4 are aligned. Then we can choose the transmitting beamformers 
[uj, u 2 ] that satisfy 



/ 



[ul,u 2 r ] CNull 



H r2 
H r 4 



Null 



d 2)3 H ri 2 
d| 3 H ri2 
d2 i4 H r2 
d2 4 H rj2 



(20) 



d 4il H r) 4 

V L ^4,lH r ,4 J J 

which is feasible since the dimension of the concatenated effective channel matrices is degraded 
to 6 x 8. The other transmitting and receiving beamformers are designed similarly so that the 
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interference streams are nulled out for each source node. Fig. |4] illustrates the process of joint 
transceiver design for interference cancellation between source nodes and the relay in the BC 
phase. 

B. Necessary Condition for d = 2N when N < 2M 

From the above subsection, we can see that the proposed scheme SAJIC achieves the DOF 
upper bound for M = 5, N = 8. In this subsection, we analyze the condition to achieve the 
DOF upper bound when N < 2M. 

In order to maximize the DOF, it can be easily seen that the number of data streams between 
each pair of communicating source nodes should be the same, i.e., cf li3 = d 3 i = d 13 , di^ = 
g?4,i — du, g?2,3 = ^3,2 — ^23, ^2,4 = ^4,2 = <^24- Note that the DOF indicates the maximum 
number of independent data streams that can be simultaneously transmitted in the network. 

In order for the relay to obtain the network coded messages, we should align the signals 
between each pair of source nodes during the MAC phase. Taking source nodes 1,3 as an 
example, we should satisfy the following conditions 



It can be seen that the dimension of the intersection space of the channels of source nodes 1 and 
3 determines the maximum number of data streams we can align. We must have the following 



H 1)r v* 3 = H 3 , r v*i ^ g r fc , 1 < k < d 13 . 



(21) 



dim (span ([gj, ...,g^ 13 ])) < dim (span (H 1>r ) n span(H 3jr )) . 



(22) 



According to the dimension theorem E51 . we obtain that 



dim (span(Hi >r ) fi span(H 3jr )) 



dim (span(Hi /r )) + dim (span(H 3)r )) — dim (span ([Hi >r H 3>r ])) 



min {M,N} + min {M, N} — min {2M, A^} 



2 min {M, N} - N 



(23) 



Combining (|22l) -d23l. we have 



dim (span ( [gj, g^ 13 ] ) ) < 2 min {M, N} - N 



(24) 
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which is equivalent to 

d 13 < 2 min {M, N} - N. (25) 
For the other pairs of communicating source nodes, we similarly have 

{d u , d 23 , d 24 ] < 2 min {M, N} - N. (26) 
Combining (f25])-(l26l), we have 

di 3 + d 14 + d 23 + d 24 <8 min {M,N}- AN. (27) 
In order to achieve the upper bound, we must satisfy the following condition 

dis + d 14 + d 23 + d 24 = N. (28) 
Based on (1271) and (1281 ), we obtain 

5iV < 8 min {M, N} 

which is equivalent to 

8M 

N< . (29) 

5 

Thus, we have obtained the necessary condition to achieve the DOF upper bound when N < 2M 
for the MIMO two-way X relay channel. 

V. ACHIEVABILITY OF THE UPPER BOUND 

In this section, we generalize SAJIC in Section IV.A. to arbitrary N, M with TV < 2M, and 
show that it achieves the DOF upper bound when iV < I . Therefore the necessary condition 
in Section IV.B to achieve the DOF upper bound is also sufficient. 

We first provide the transmission scheme for the case of N = ^,VM = 5k, k 6 Z + and 
show that (d 1)3j d 1>4 , d 2 , 3 , d 2A , d 3A , d 3)2 , d 4>1 , d 4>2 ) = (f , f , f , f , f , f , f , f ) is achieved by this 
scheme. 

During the MAC phase, the ith source node sends message Wij to the jth source node using ^ 



independently encoded streams along beamforming vectors Vj j 
the transmitted signal from source node 1 is 

x l = » 1,3^1,3 + Vi^Si^ 



JV 



V 1 ■ V 4 . 



. For instance, 



N_ 

4 



E «A + (30) 



k=l 
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where s x 3 and s x 4 are the jX 1 encoded symbol vectors for W\$ and W\^, respectively. The 
transmitted signals from other source nodes are in a similar form. In order for the relay to obtain 
the network coded messages W 1)3 © W 3j i, Wi,4 © W^i, © ^3,2 and W 2 ,4 © W it2 , we should 
carefully choose the beamforming vectors to satisfy the following signal alignment conditions 



Hi, r vf i3 


TT - r fc A 


gr ' 






= H^V*^ = 


gr , 




H2,rV2 )3 


= H 3jr V3 2 = 


f+fe 




H2,r v 2,4 = 


- H 4jr V4 2 = £ 


5»* ; 


1 < k < 



N 

T 



(3D 



where 



are iV transmitting vectors seen by the relay. The above conditions imply that 



span 
span 
span 
span 



JV 



- + 1 



AN 

lr 4 



C span (H ljr ) fl span (H 3iT .) 
C span (Hi )r ) fl span (H 4ir ) 
C span (H 2 , r ) fl span (H 3)J .) 
C span (H 2 r ) fl span (H 4jr ) . 



jv 



(32) 



Since all the entries of the channel matrices are i.i.d. zero-mean complex Gaussian random 
variables, there exists a (2M — N = ^) -dimensional intersection subspace constituted by the 
column space of channel matrices for each pair of communicating source nodes with probability 
1. Then we can always choose ^ linearly independent transmitting vectors {g^ } for each source 
node pair. As a result, the received signal in © is rewritten as follows 



n, 



(33) 



JV 
4 



where the N x N matrix G r = [g£, g^] , and the N x 1 vector s r = [s\ 3 + s 3 1; s^ 3 + 

S 3*l' 5 1,4 S 4,l' S l*4 "I" "^4*1 ? ,S 2,3~'~ S 3,2> "•) S 2^3 S 3,2' S 2,4 S 4,2' "•> S 2*4 ■ AlsO since the 

entries of all channel matrices are independently Gaussian, the probability that a basis vector 
in the intersection space of one pair of source nodes' channel matrices lies in the intersection 
space of another pair is zero. Thus G r is full-rank with probability 1, which guarantees the 
decodability of s r at the relay. The four network coded messages W\ 3 = Wis © W 3 i,Wu = 
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Wi t 4 © W 4) i, W23 = W 2 ,3 © W 3t2 and W 2 4 = W 2) 4 © are then obtained by applying the 
mapping principle of physical layer network coding [2] to each entry of s r . 

For the BC phase, the relay broadcasts the network coded messages W13, Wu, W23 an d 
W24 to all source nodes using encoded symbols q r = [ql, q^] T along the beamforming 

£L — + 1 — — +1 — 

vectors U r = [u*, uf ]. More specifically, [q^, g r 4 ] T , [g r 4 , q r 2 ] T , [q r 2 , q r 4 } T and 
[g r 4 , g^] T are the ^ x 1 encoded symbol vectors for FT 13 , W14, W23 and W24, respectively. 
Then the transmitted signal at the relay in © is rewritten as 

N 

x r = ^u r fe <^. (34) 

k=l 

The received signal at source node 1 is given by 



yi = H ri ix r + ri! 

/ f ¥ n x 

= D x H r)1 (J2 u r<lr+ E + ^ u^ r fc +n a (35) 

V fc=i^ / fe=f+i ^=^+1 / 

signal interference 
where the first term in the bracket represents the combination of the desired network-coded 

messages W13 and W14, while the remaining two terms are the unwanted interference W23 and 

W24. The received signals at the other source nodes are written in a similar way. We can see 

that each source node suffers from two sources of interference, i.e., the signals intended for the 

other two pairs. 

Next, we jointly design the transceivers for the source nodes and the relay for interference 
cancellation. For the source nodes, we design their receiving matrices to align the receiving 
effective channels for each pair of source nodes as follows 



dl,3Hf-,l 


= d^iHr^ 


A_ 


W ^3> 


di j4 H rjl 


= d 41 H r) 4 


_A 


W 14> 




= d3 j2 H ri 3 


A_ 


W23, 


d2 )4 H rj 2 


= d 42 H r . ) 4 


_A 


W 24, 



1 < k < f (36) 



where wf is a 1 x N effective channel vector. The above conditions are equivalent to the 



14 



following 



span f 


1 T 
W 13 >■ 


N T- 


span f 


1 T 
w 14 ,. 


nT~ 
■»Wi4 


span ( 


i r 

w 23 >■ 


jvT" 
•> w 23 


span f 


i r 
w 24 ,• 


jvT" 

-,w 2 4 4 



C span (HjTjJ D span (H^ 3 ) 

C span (H^ x ) fl span (H^ 4 ) 

C span (H£ 2 ) n span (H£ 3 ) 

C span (H^ 2 ) n span (H^ 4 ) . 



(37) 



Similar to the MAC phase, there always exists a (2M — N) = ^-dimensional intersection 
subspace constituted by the row space of the channel matrices for each source node pair with 
probability 1. Then we can choose iV linearly independent vectors to satisfy the above conditions. 

For the beamforming vectors at the relay, we can choose them to lie in the intersection subspace 
of each source node pair's effective channels' null space as follows 

H r 4 



span 



u. 



r , . . . , U r 



C Null 



span 



span 



span 



u, 



u, 



u, 



L El 

, ...,u r 2 



-+1 



u r 4 



.N 



C Null 



C Null 



C Null 



H ri2 

H r ,3 

H r> i 
H n i 

H r ,3 



(38) 



We show that there exists a ^-dimensional null space for the concatenated effective channel 
matrix of each pair of communicating source nodes with probability 1. Taking source nodes 

T 



2 and 4 as an example, the dimension of 
their receiving effective channels in d36l) . 



I I / ttT 

n r,2) il r-,4 

H~ T tt7 
t.2-1 n r,4 



is x N. Since we have aligned 
has ^ repeated rows and its rank is 



min{iV — ^j, iV} = ™. Therefore, the dimension of its null space is iV — ^ = ^. For the the 
other source node pairs, we can similarly get the result. Also due to the fact that the null space 
of the concatenated effective channel matrix for each pair of source nodes has no intersection 
with those of the other pairs, all iV beamforming vectors at the relay are linearly independent 
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with probability 1. Thus the received signals at source node 1 can be rewritten as 



JV AT 

yi = H r /]Tu r \ fc + Yl <Qr)+^- (39) 

^ k=i fe=f +i 



for ^ for Wli 

Thus, there is no interference for source node 1 and it can decode these useful signals. Then 
by using its own messages, source node 1 can obtain the messages from source nodes 3,4 as 
follows 

W 3) i = W h3 © W l3 , W AA = W 1A © W u . (40) 

In the same manner, the other source nodes can also obtain the messages intended for themselves. 
Therefore, a total of 2N DOF is achieved by using the proposed scheme on MIMO two-way X 
relay channel. 

A. N < [^f\ or M ^ 5k 

For the other cases that M ^ 5k or N < [^J , we can choose the DOF for each pair for 
different values of N as below 

N = 4k: 

4,3 = 4,1 = X> 4,4 = 4,1 = f 
4,3 = 4,2 = X' 4,4 = 4,2 = ^ 

N = 4k + 1 : 

4,3 = 4,i = LfJ > rf M = d ^ = LfJ 

4,3 = 4,2 = LfJ . 4,4 = 4,2 = Lf J + 1 

N = Ak + 2 : 

4,3 = 4,i = LfJ . 4,4 = 4,i = Lf J + 1 

4,3 = 4,2 = Lf J + 1, 4,4 = 4,2 = LfJ 

N = 4k + 3 : 

4,3 = 4,1 = LfJ . 4,4 = 4,i = Lf J + 1 

4,3 = 4,2 = LfJ + 1, 4,4 = 4,2 = Lf J + 1 



16 



We can similarly apply the previous transmission scheme to achieve the upper DOF bound 2N 
and the process is briefly described as follows. 

For the MAC phase, we show that the signals for each pair of source nodes can be aligned 
at the relay: since 



N < 



8M 



< 



8M 



we can have 



2M -N> 



when N = Ak 



(41) 



LfJ +1, when N ^ Ak. 
For the BC phase, the receiving alignment is also feasible just as the MAC phase. For the 
transmitting beamforming design at the relay, we show that relay can always choose linearly 
independent beamforming vectors for each part of the signals. Without loss of generality, we 
take [uj, ...,u^? 13 ] as an example, which should satisfy the following condition 

H ri4 

Since 



span([ U ;,...,< 3 ]) C Null 



(42) 



H 



r,2 



w 23 




wi 4 


w 23 
W 24 


j H ri4 — 


w 14 

W 24 




. W 24 




. W 24 



(43) 



we have 




= d 13 = dm([uj,...,u* s ]). 
Finally, we summarize the algorithm for SAJIC in the following chart 



(44) 



Outline of SAJIC 
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• Step 1. In the MAC phase, each source node designs its beamforming vectors {vf-} 
according to (13D) so that the two desired signals for network coding are aligned at the relay 
node. 

• Step 2. By applying the mapping principle of physical layer network coding, the relay 
then decodes its received signals to obtain the network coded messages {Wi?}. 

• Step 3. In the BC phase, the source nodes and the relay jointly design their transceivers. 
More specifically, all the source nodes design their receiving filter matrices {D^} according 
to d36l) : the relay designs its transmitting beamforming vectors {u^} according to (1381) to 
cancel the interference for each source node. 

• Step 4. Each source node decodes its received signals to obtain the network coded messages 
{Wij} intended for itself. Using its side-information, each node finally acquires its desired 
messages {H^j}. 

Remark 1: We can see that the DOF for each source node may not be the same when N is 
not a multiple of 4. However, we can apply time slot extension here to let every source node 
achieve the same DOF y lfi~5l . For instance, in the case M = 5, N = 7, in the first time slot, 
we let 



In the second time slot, we let 



In time slot 3, we let 



And in time slot 4, we let 



4,3 


= 4,i 


= 1, 


4,4 — 


4,1 


= 2. 


d 2 ,3 


= 4,2 


= 2, 


4,4 — 


4,2 


= 2 


4,3 


= 4,i 


= 2, 


4,4 = 


4,1 


= 1, 


4,3 


= 4,2 


= 2, 


4,4 — 


4,2 


= 2 


4,3 


= 4,1 


= 2, 


4,4 = 


4,1 


= 2. 


4,3 


= 4,2 


= 1, 


4,4 — 


4,2 


= 2 


4,3 


= 4,1 


= 2, 


4,4 = 


4,i 


= 2. 


4,3 


= 4,2 


= 2, 


4,4 — 


4,2 


= 1 
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Thus, each source node can achieve the same DOF of |. 

In Section IV.B we have shown that N < \ is a necessary condition to achieve the DOF 
upper bound. And in this section we have shown that this condition is also sufficient. Hence we 
have the following main result of this paper. 

Theorem 2: When iV < 2M, the necessary and sufficient condition for SAJIC to achieve the 
DOF upper bound 2N in the MEMO two-way X relay channel is iV < |_^J . 

Remark 2: If channel reciprocity holds, we can see that the beamforming vectors of the source 
nodes {v^} in the MAC phase are the same with their receiving filter vectors {d^ T } in the 
BC phase. 

Remark 3: When [_^rJ < N < 2M, SAJIC still works but no longer achieves the upper 
bound 2N. When iV > 2M, the SAJIC will not work since there exists no intersection subspace 
between each pair of source nodes' channel matrices and therefore signal alignment is no longer 
feasible. 

VI. Connection with the Transmission Method in [|23Tl 

In the previous sections, we have shown that using our SAJIC, the cut-set outer bound for 
the DOF can be achieved when N < |_^J . Specifically, we align the signals for each pair of 
source nodes in the MAC phase and apply joint transceiver design for interference cancellation 
in the BC phase. In this section, we will show that if we do not consider the joint transceiver 
design but directly apply interference nulling beamforming at the relay in the BC phase, our 
proposed scheme will reduce to a generalized version of the transmission method in [|23l . 

For the reduced or simplified transmission scheme which does not apply the joint transceiver 
design in the BC phase, we consider as an example the case N = ^,VM = 3k, k G 
Z+, (d lj3 , d 1A , d 2i3 , d 2A , d 3)lj 4,2, d 4 ,i, M = (f , f , f , f , f , f , f , f )• In the MAC phase, we 
similarly apply signal alignment as in (|3T| ). In the BC phase, since we directly apply interference 
nulling at the relay, the relay will cancel one part of interference for each source node, leaving 
the remaining part of interference to be cancelled by the source node itself. More specifically, 
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the relay can choose its beamformers as 



span 
span 
span 
span 

For each channel matrix EL,-, i 



u 



u 



-+i 



-+i 



u 



u 



3N -\ 

4 



C Null(H P)4 ) 
C Null(H r , 2 ) 
C Null(H nl ) 

CNull(EL j3 ). (40) 
1, ...,4, there exists a (N — M = ^) -dimensional null space 



u r 4 , ...,u r 



N 



with probability 1. Then for each network coded message, the relay can choose ^- linearly 
independent vectors. Also it can be easily seen that the iV beamforming vectors are linearly 
independent with probability 1. Thus the received signals at source node 1 can be rewritten as 



JV 

yi = H ril I U r Qr + 

^ k=l 



signal 




+ 11!. 



(41) 



interference 



Note that source node 1 has M = ^f- antennas and the dimension of the useful signal is —. So it 



has exactly ( 



37V 
4 



N 
2 



) free dimensions for the interference signal whose dimension is also 



; j. Thus source node 1 can cancel the other part of the interference by itself. More specifically, 



N n t 

source node 1 can choose its receiving matrix Di 6 C~ as follows 



span (Df ) C Null 



:',.V 
4 



H r>1 u r 4 , H rjl u 



N 



(42) 



Since 



:'.,V 
4 



,H u 



N 



isanMx j matrix, the dimension of its left null space is M — ^ 



Y- Source node 1 can choose y linearly independent receiving filter row vectors, and hence 
is full-rank with probability 1. Then the received signals for source node 1 is 



yi = Di 



4 



2 




(43) 



for Wi 3 



for w 14 



There is no interference and source node 1 achieves the DOF of y. The other source nodes 
operate in the same manner. Therefore, the reduced transmission scheme also achieves the total 
DOF of 2N. 

Remark 3: The method given in E3l corresponds to the reduced transmission scheme for 

N = 4,M = 3. 



20 



Next, we show that the reduced transmission scheme requires a stricter condition to achieve 
the DOF upper bound, i.e. N < [^\ . 

In the MAC phase, the condition that the reduced scheme needs to satisfy is 

8M 



N < 



(44) 



5 

which is the same as the original proposed transmission scheme. Extra conditions are needed in 
the BC phase for the reduced scheme. 

Lemma 1: For the reduced transmission scheme in the BC phase, for each source node, the 
dimension of the interference that needs to be canceled by the relay iJ3 N — M. 

Proof: Without loss of generality, we take source node 1 as an example. For source node 
1, the dimension of useful signals is diz + du\ while the dimension of interference is d 23 + ^24- 
Since it has M antennas, the interference dimension that it can cancel by itself is M—(di3+du). 
Then the interference dimension that needs to be nulled at the relay is 

d 23 + d 2 4 — [M — (d 13 + d u )} 
= d 13 + du + d 2 3 + d 2 4 — M 
= N-M. 

For the other source nodes, we can similarly obtain the result and the lemma follows. ■ 
For source node 1, suppose d\ z out of d 23 interference streams and d\± out of d 2i interference 
streams are nulled out at the relay. Then according to Lemma 1, we have 

4s + d\ A = N - M. (42) 

Since the dimension of the null space of H r A is also iV — M, the relay can choose beamformers 
which lies in its channel matrix's null space to cancel these interference streams. For source 
node 2, we can similarly have that 

d\ 3 + d{ 4 = N-M. (43) 
We now consider source node 3 and source node 4. 

'if N < M, the relay does not need to cancel the interference and each source node can null all the interference it suffers 
by itself. 



21 



Lemma 2: During the BC phase, the null space of each source node's channel matrix has no 
intersection with that of the other source nodes, i.e., 

Null(H ni ) n Null(H rj ) = 0, Vz ^ j. (44) 

Proof: We first consider the channel matrices for source nodes 1 and 3 and have that 



Null(H r>1 ) n Null(H r . 3 ) = Null _ 

l T 



(45) 



Since [H^HjTg] is a 2M x A" matrix and A" < 2M by assumption, the dimension of its 
null space is always zero. For the other pair of source nodes, the same argument holds and the 
lemma follows. ■ 
As for source node 3, the interference signals consist of du + d 24: data streams. From the 
previous discussion, we know that d\ 4 out of d\± data streams lie in the null space of source 
node 2's channel matrix; and d\ A out of d 2i data streams lie in the null space of source node l's 
channel matrix. According to Lemma 2, these interference signals cannot lie in the null space 
of H r 3 and therefore cannot be nulled out by the relay. Thus, source node 3 must cancel them 
by itself, which implies that 

d\ A + 4 4 < M - (di 3 + d 23 ) ■ (46) 
Similarly for source node 4, we should have 

d\ 3 + d\ 3 <M- (d 14 + d 24 ) • (47) 
Combining (l46l) and (l47l) . we have 

d{ 3 + d\ A + 4 3 + 4 4 <2M - (d 13 + d u + d 23 + d 2i ) . (48) 



(49) 



Plugging d28j), d42j) and (g3j) into (|48|). we then obtain 

AM 



N < 



3 

which is the condition the reduced scheme should satisfy in the BC phase. Now combining (1441) 
and (|49l , the necessary condition for the reduced scheme to achieve the upper bound becomes 
simply (|49l . 

From the above analysis, it can be seen that using the reduced transmission scheme, we 
can achieve the DOF upper bound 2Af in a range < A^ < I . By further applying joint 



22 



transceiver design for interference cancellation in the BC phase, our proposed scheme can achieve 
the upper bound 2iV in a wider range < N < [^f\ ■ 



VII. Simulation Results 

In this section, we provide numerical results to show the ergodic sum rate performance for 
the proposed transmission scheme. Then, we will demonstrate that the proposed scheme exactly 
attains the upper bound on the DOF derived in Section V. The channel is assumed as the 
normalized Rayleigh fading channel, i.e., the elements of each channel vector are independent 
and identically distributed circularly symmetric zero-mean complex Gaussian random variables 
with unit variance. The numerical results are illustrated with respect to the ratio of the total 
transmitted signal power to the noise variance at each receive antenna in decibels (SNR = P). 
Each result is averaged over 10000 independent channel realizations. 

We now explain how we compute the sum rate for the MIMO two-way X relay channel 
when applying the SAJIC. In the MAC phase, assuming that the zero forcing detector W r = 

-i T 

rp rp -t a 

is applied by the relay, the achievable rate for network coded message Wis is 



calculated as 

#13 



log 



det I + 



1 T 



H 



w 



di3 



H 



In the BC phase, the achievable rate for Wis at source node 1 and 3 is given by 



R' 13 = lo£ 



det 1 + 



4 T 



T 



w 13 ,».,W 13 



u r uf 



H 



Then we have 



#1,3 = #3,i = min {R 13 , R' 13 } . 



(50) 



(51) 



(52) 



For the rates of other pairs of communicating source nodes, we can compute them in a similar 
way. Thus, we can obtain the achievable sum rate for the whole network when SAJIC is applied. 

In Fig. [5l we plot the sum rate performance of the proposed scheme according to various 
antenna configurations. We can see that, as analyzed in Section V, the proposed scheme indeed 
achieves the upper bound on the DOF. In specific, we can always observe a sum-rate increase 
of 27V bps/Hz for every 3 dB increase in SNR. For instance, when M — 5, N — 8, the curve 
has a slope of 2 N = 16. In Fig. [6l we plot the sum rate performance of the network when the 
number of antennas at each source M is fixed. It can be seen that as the number of antennas at 
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the relay N increases, the total DOF also increase, which shows that the relay antenna number 
is a bottleneck of the network when N < 2M. 

VIII. Conclusion 

This paper considered the total DOF for the MIMO two-way X relay channel. We analyzed 
the upper bound on the DOF for such a network. Then by exploiting physical layer network 
coding and joint interference cancellation, we proposed SAJIC and showed that SAJIC can 
achieve the upper bound if and only if N < . Also, we generalized the scheme in ||23l 
and derive a necessary condition for it to achieve the upper bound. Besides, we analyzed the 
relationship between these two schemes and showed that our proposed SAJIC can be reduced 
to the generalized version of the scheme in E3l after some simplification. 

The achievability of the upper bound on the DOF for the considered network in the case of 
L^NpJ < N < 2M and the case of N > 2M remains open for further investigation. 
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Fig. 1. MIMO two-way X relay channel. 




Fig. 2. One direction of network information flow for the MIMO two-way X relay channel. 
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Fig. 3. Signal alignment for network coding during the MAC phase. 
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Fig. 4. Joint interference cancellation during the BC phase. 




Fig. 5. The ergodic sum rate for the MIMO two-way X relay channel under different network architectures. 
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Fig. 6. The ergodic sum rate for the MIMO two-way X relay channel when M is fixed. 



